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Abstract 

The analytic solution of Teukolsky equation in Kerr-de Sitter and Kerr-Newman- 
de Sitter geometries is presented and the properties of the solution are examined. In 
particular, we show that our solution satisfies the Teukolsky-Starobinsky identities 
explicitly and fix the relative normalization between solutions with the spin weight 
s and —s. 
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1 Introduction 



In our series of papers, Mario, Suzuki and Takasugi [1, 2] have constructed the ana- 
lytic solution of the perturbation equation of massless fields in Kerr geometries which is 
commonly called Teukolsky equation [3]. The solution enabled us to investigate various 
properties of black hole analytically [4], the scattering problem of particle emitted in black 
hole, the analytical expression for the absorption rate by Kerr black hole which was one 
of the main concern by Chandrasekhar [5], and the Teukolsky-Starobinsky identities. 

Our solution is expressed in the form of series of hypergeometric functions and Coulomb 
wave functions. The coefficients of these series are determined by solving three term re- 
currence relations. Two series have different convergence regions and they are matched 
in the region where both series converge. This method turns out to be quite powerful 
in a practical calculation and is sucessfully applied to examine gravitational waves from 
a binary in a post-Newtonian expansion; to construct the template of the gravitational 
wave emitted from a particle moving around a Kerr black hole [6] and the rate of the 
gravitational wave absorbed by the black hole [7]. 

In our previous paper [8], Suzuki, Takasugi and Umetsu have extended this method 
to solve the perturbation equations of massless fields (the Teukolsky equation) in Kerr-de 
Sitter and Kerr-Newman-de Sitter geometries. We found the transformations such that 
both the angular and the radial equations are reduced to the Heun's equation [9]. The 
solution of Heun's equation is expressed in the form of series of hypergeometric functions 
and its coefficients are determined by three term recurrence relations, similarly to the 
Kerr geometry case. It should be noted that electromagnetic field and gravitational fields 
couple each other in Kerr-Newman-de Sitter geometries and do not follow the equation 
we considered there, although the fields follow the equation in Kerr-de Sitter geometries. 

The solution of the radial Teukolsky equation given in Ref. [8] is convergent for r < r' + 
with r' + being the de Sitter horizon. In this paper, we give the solution valid in the entire 
physical region. This is achieved by constructing another solution which is convergent for 
r + < r and then by matching this solution and the previous soluiton in the reigion where 
both solutions are convergent. Then, we examine porperties of the solution in detail. 
In particular, we show analytically that our solution satisfies the Teukolsky-Statrobinsky 
(T-S) identities [10], if we properly take the relative normalizaition between the solution 
with the spin weight s and the one with — s. 

It is quite surprising that the analytic solution is obtained for Kerr-de Sitter and Kerr- 
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Newman-de Sitter geometies, similarly to the Kerr geometry case. This might be the 
reflection of an underlying symmetry which the gravitational theory contains. Although 
the use of this solution is less clear in comparison with that in the Kerr geometries, we 
believe that our finding is important not only because the technique we used can be 
exteded to other phenomena, but also because the solution may become relevant to the 
early stage of universe where the cosmological constant play important role. 

In Sec. 2, we review the transformation of the radial equation to Heun's equation, 
because we change some definitions of parameters from the previous paper [8]. In Sec. 3, 
the solution is given in the form of series of hypergeometric functions. In particular, 
we present the solution which is convergent around the de Sitter horizon. We discuss the 
convergence region of the solution. Then, we match two solutions in the region where both 
solutions converge. In Sec. 4, we explicitly show that our solution satisfies the Teukolsky- 
Starobinsky identities and the relative normalization between the solutions with the spin 
weight s and —s. We found some identities involving coefficients of series. Summary and 
discussions are given in Sec. 5. 

2 Teukolsky equation for the Kerr (-Newman) -de Sit- 
ter geometry and Heun's equation 

We consider the Teukolsky equations for the Kerr-Newman-de Sitter geometries. In the 
Boyer-Lindquist coordinates the Kerr-Newman-de Sitter metric has the form, 

,2 2 ( dr 2 d9 2 \ A e sin 2 6> r , , 9 9x , l9 

ds = - p {x + ^)-jrTaf? [adt - {r+a)M 

A ' (dt-asm 2 9d<p) 2 , (2.1) 



1 + a) v 



where 



A r = (r 2 + a 2 ) ( 1 - %r 2 ] - 2Mr + Q 2 = -%{r - r+)(r - r_)(r - r')(r - r'_), 
V a 2 J a 2 

A 9 = 1 + a cos 0, a = — — , 

p = r + ia cos 9, p 2 = pp*. 

(2.2) 

Here A is the cosmological constant, M is the mass of the black hole, aM its angular 
momentum and Q its charge. 
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We assume that the time and the azimuthal angle dependence of the field has the form 
e -i(wt-m4) _ Then, the equation for the radial part with spin s and charge e is given by 

2a 



A r V{D\ + 2(1 + a) (2s - l)iur - — (s - l)(2s - 1) 

a 

-2(1 + a)eQKr + iseQrd r A r + e 2 Q 2 r 2 
H 2iseQ — s(l — a) — A s 



R s = 0, 



(2.3) 



where 



t(l + a)K d r A 
V n = d r h n 



A r 



A r ' 



VI = dr + ^ 1 + ^ K + n d ^ 



(2.4) 



with K = u(r- 2 + a 2 ) — am. In our previous paper [8], we showed that solutions are 
characterized by the characteristic exponent (shifted angular momentum) v. 

In this paper, we use slightly different definitions of the separation constant \ s and the 
characteristic exponents v and — v — u, where u is the parameter appearing in the Heun's 
equation, is given explicitly in Eq.(2.17) and is different from the angular frequency in 
Eq.(2.3). We redefined \ s and v as 



A, 



A s + s{\ - a) , 



(2.5) 



With these modifications, as we showed in Ref.[8], A s becomes an even function of spin 
variable s, 



A s — X-s , 



(2.6) 



and if v is a characteristic exponent, then —v — \ becomes the one. That is, both 



v and — v — 1 . 



(2.7) 



become characteristic exponents. 

With a new A s , the radial Teukolsky equation is explicitly written by 

eQr \ 2 . . / eQr \ dA r 



. _, d . , , , d 1 

A — A s+1 1 

r dr r dr A, 



r 

2a 

2 



dr 



+ Ais(l + a)ujr -(s + l)(2s + l)r 2 + s(l - a) - 2iseQ - \ s }R S = 0. 



a 



(2.1 



This equation has five regular singularities at r±,r' ± and oo which are assigned such that 
r± -> M ± VM 2 - a 2 - Q 2 and r+ -> ±-^= in the limit a -> (A -> 0). 

Next, we use a variable x rather than z used in our previous paper [8], which are 
defined by 

x = 1-^ = ) ^7 rf 2.9 

This transformation from r to x maps the inner horizon r_, the outer horizon r+, the de 
Sitter horizon r' + , r'_ and oo to 0, 1, and Xoo, respectively, 

_ (r_ - rL) « - r+) 
(r_ - r+) (r+ - r_) 

Xoo = l-z OQ = ) =(. (2.10) 

In the followings, we omit the subscript s for the radial wave function. 
To proceed further, we define the following parameters, 

A i± = ±{- s ±(2a t + s)} (i = l,2,3,4), (2.11) 

Here a's are pure imaginary numbers defined by 



ai 



.a 2 (l + a) (uj{r\ + a 2 ) 



eQr+ 



a 2 = i 



a (r' + — r+)(r'_ — r + )(r_ — r+) 
a 2 (l + a) (u{rl + a 2 ) - am - 



a 3 



a (r' + — r_)(r!_ — r_)(r + — r_) 

,a 2 (l + q) [ u ( r + + a )- a ^-T+^j 
a (r_ -rV)(r'_ -r^)(r+ - /+) ' 



.a 2 (l + a) (4' 2 + « 2 )-™-T 
and they satisfy the relation 



°4 = * 7 TTT^ TT ( T\ ( 2 - 12 ) 

a (r_ — r_){r + — r_)(r + — r_) 



CL\ + 02 + Gt3 + 04 = 0. (2-13) 

Now, following the transformation given in Ref.[8], we can factor out the singularity 
at x^ by the transformation as 

K„ ;{ o, 1};s (-) = (-^(l-^(^f) (^) W*), (2.14) 
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where we choose A\ = Ai_ for the solution to satisfy the incoming boundary condition 
at the outer horizon, and the other Ai(i = 2,3) takes either A i+ or A^. Then, we find 
that /j'g X | (x) satisfies the following Heun's equation 



dx 2 



+ 



7 5 e 
- + 7 + 



X X 1 X Xf 



d a + a_x + v 

1 - 

dx x(x — l)(x — x r 



} /{W*) = > ( 2 - 15 ) 



where 



7 = 2A X + s + 1, 5 = 2A 2 + s + 1, e = 2A 3 + s + 1 , 

a± = A x + A 2 + A 3 + A 4T + 2s + 1 . 

Parameters 7, 5, e and er± satisfy the following relation, 

7 + 5 — 1 — a + + a- — e = a; , 



(2.16) 



(2.17) 



which is required for Eq.(2.15) to be a Heun's equation. It should be noted that the pa- 
rameter u defined above is a commonly used notation for Heun's equation and is different 
from the angular frequency. The remaining parameter v is given by 

2a 4 (l + a) 2 (r_ - <) 2 (r_ - r'_)\r + - r'_){r' + - r'_) 



v = 



a 2 V 



r + — r_ 



l — LU 2 r+(r + r_ — 2r_r' + + r + r' + ) + 2au(auJ — m)r + {r_r' + — r 2 + ) 

2 / \ 2 / / \ r 2 / 2 / f \ 

—a (aw — to) (2r + — r_ — r + ) + :j— ^ — [ wr + (r + r_ + r + — 3r_r + + r + r + ) 
— a(acj — m)(r + r_ — 3r 2 + + r_r^ + r+r^) ] + r +( — r + + r - r ' + ) j 

2isa 2 (l + a) 



cj(r + r'_ + a 2 ) — am — -£g- r++r - 



-a 2 



(r+ — r_)(r+ — r'_)(r + — rl) 



- (s + 1) ^(2s + 1) 

- 2A 1 (x r A 2 + A 3 ) + 



2r' 



(r+ — r_)(r' + — r'_) 



+ [{l + x r )A 1 + x r A 2 + A 3 ] 



a(r + — r_)(r' + — r'_) 
where T> is the discriminant of A r = 0, 



— \-\ s + s(l - a)} , 



(2.18) 



V 



(r+ - r_) 2 (r + - r' + ) 2 (r + - r'_ ) 2 (r_ - r' + ) 2 (r_ - r'_) 2 (r' + - r'_f 



16a 



10 



a- 



{l-af M 2 - (1 - a)(a 2 + Q 2 



16a 2 



2\3 



+ - 



ft 



(a 2 + g 2 ) 



-27M 4 + 36(1 - a)M 2 (a 2 + Q 2 ) - 8(1 - a) 2 (a 2 + Q 2 ^ 2 



(2.19) 
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A parameter u in the definition of v is an angular frequency. 

Modifications in Eq.(2.5) and the use of the variable x rather than z cause changes of 
parameters in the Heun's equation from those in Ref. [8] . The changes of parameters are 
shown below: 

B 1 -> A 2 , B 2 -> A u B 3 -> A3, 
(a, (5) — > (cr+, (7_), — g — > i>, t> — > — i> — o" + cr_, a# — > x r . (2.20) 

Parameters in the left-hand side are those used in Ref. [8] and those in the right-hand 
side are those used in this paper. The v in Eq.(2.18) is derived from the corresponding 
parameter v of Eq.(3.10) in Ref. [8] by exchanging r + and r_ and changing of \ s in Eq.(2.5). 
There was a miss-print in v in our previous paper [8] where — \ s — 2iseQ should be read 
-Xs- 

3 Solutions of Teukolsky equation 

In this section, we derive the solutions which are convergent in the region including the 
outer and the inner horizon and also the solutions which are convergent in the region 
including the de Sitter horizon and 00. Then, by matching these solutions in the region 
where both solutions are convergent, we obtain the solutions valid in the entire region. 

(A) Solutions convergent around outer horizon 

There are two independent solutions; one satisfies the incoming boundary condition 
at the outer horizon and the other does the outgoing boundary condition. They are 
expressed as series of hypergeometric functions and these series converges in the region, 
r < r' + with r' + being the de Sitter horizon. 

(a-1) The solution which satisfies the incoming boundary condition at the outer horizon 
The solution is given by taking into account of changes in Eq.(2.5) and adopting the 
variable x as 

(rp rp \ -^4.3 / <T™ y» \ 2s - ) - ! 

f{o,iy,,( x ) = £ < F (-^- z/ + |-^ n + z/ + | + ^7;^) , (3.1) 

n=— 00 ^ ' 

where v is the characteristic exponent (the shifted angular momentum). Coefficients a v n 
are determined by solving the following three term recurrence relation with the initial 
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condition = 1, 



«X+i + « + «-i = 0, (3.2) 

where 

(n + z/-f + |)(n + z/- ( 7 + + f + |)(n + z/- ( 7_ + f + |)(n + z/ + ^- f + |) 

2(n + u+ l)(2n + 2z/ + 3) 

1 CiJ 2 — 1 

+^K7 - 5) + 5(1 - w) + 2<j + ct_] + — - — x r + 

(w + i/ + o-+ - f - ^)(w + i/ + o-_ - f - i)(w + i/ + 7 - f - i)(w + i/ + f - ^) 
7n - 2(n + i/)(2n + 2i/-l) 

(3.3) 

Determination of the characteristic exponent 
The characteristic exponent v is determined such that the series of hypergeometric 
functions converges. Let us define the continued fractions 



Rniy) = Lniy) = (3.4) 

which satisfy 



a n-l "n+1 



T or 

Now we can evaluate the coefficients a v n by using either series Rniy) or L n [y) with an ap- 
propriate initial data. The convergence of the series requires the following transcendental 
equation, 

R n (is)L n ^(v) = 1. (3.6) 

This is the equation to determine the characteristic exponent v. By using aC^ -1 = 7^ 
and Pz^ 1 — Pni we can P rove that if z/ is a solution of the transcendental equation, then 
— v — 1 is a solution, similarly to the proof given in Ref.fl] for the Kerr geometry case. 
We can also prove that 

air 1 = <C (3-7) 
by using the recurrence relation and by taking the initial condition a v § = a^" -1 . 
(a-3) The solution which satisfies the outgoing boundary condition at the outer horizon 
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This solution is simply given by 



Rout;{0,l};s ~ K^r * Rin;{0,l};-s) ) 



(3.J 



which is proved in Appendix. The convergence region of this solution is the same as that 

for R in;{0,l};s- 

(a- 4) Solutions specified by characteristic exponents 

By using the relation in Eq.(3.7), we can show that R^n-{oiy s ( x ) = Rin;{o,i};s( x )- This 
means that R" n .{ Q iy s ( x ) is expressed by a sum of two independent solutions with charac- 
teristic exponents, v and —v — 1. We can explicitly show with z = 1 — x that 



Rin;{0,l};s( x ) ~ R{0,l};s( z ) + R{0,l};s( Z ) 



(3.9) 



where 



iW*) = zA ^ z - l ^ Al { l -v) A \ l ~t 



2s+l 



r( 7 ) 



oo 



T(2n + 2z/+ 1) 



1 



,n+u-% + -. 



xF 



uj 1 

-n-i/+---,-n-i/ + 7 



^ 1 „ 1 

2 _ 5 ; ; ; 



(3.10) 



(a-5j T/ie convergence region of the series for R^ n .^ 01 y s (x), etc. 

Firstly, we consider the convergence region for Ri n; {o t iy >s (x) ■ We find that 



lim 

n^oo 



F» 



lim 

n^— oo 



F u 

n 



max{|e^|, |e ^|} > 1 , 



(3.11) 



where we used the abbreviated expression F% for the hypergeometric function appeared 
in Eq.(3.1) as 



F u = F 



LU 1 LU 1 

-n-v+---,n + u + - + -; T ,x 



(3.12) 



and 



= 1 - 2x ± ^(1 - 2x) 2 - 1 . 



(3.13) 



On the other hand, the ratio of coefficients converge either e^ r > 1 or e ?r < 1. Taking 
account for Eq.(3.11), we have to require their limits as follows; 



lim 



hm — — 



ra-»-oo a 



(3.14) 



n+l 



where 

S r = 1 - 2x r + - 2x r ) 2 - 1 > 1 (x r <0). (3.15) 
In order for the series of hypergeometric functions converges, it is required that 

max{|e ? |, |e^|} < e ir . (3.16) 

Now we parametrize x as 1 — 2x = (t + t)/2 where |t| > 1. This transformation maps 
a circle in the complex t plane, \t\ = C (C:constant) to an ellipse in the complex x plane. 
The condition of the convergence in Eq.(3.16) is expressed by \t\ < e ?r , i.e., inside the 
circle. In the complex x plane, the convergence region is inside the ellipse with foci and 
1, where the major axis is along the real axis extending from x r to 1 — x r . 

The convergence region for Ro U t : {o,i};s( x ) is the same as Ri n -io : n ;s (x) as seen from 
Eq.(3.8). The convergence region for R{ 0j n ;s (z) is also the same as Ri n {o,i} ;s ( x ) i because 
R{o,i};s( z ) and R7oiy 8 (z) are expressed by linear combinations of R" n .{ Q iy s and R v out .^ 1 j.. s . 

In summary, the convergence region of series of all solutions is inside of the ellipse 
with foci and 1, where the major axis is along the real axis extending from x r to 1 — x r . 
If we confine x to the physical region, the convergence region is x > x r , i.e., r < r' + , with 
r' + being the de Sitter horizon. 

(B) Solutions convergent around the de Sitter horizon 

Here, we use the variable z rather than x and construct solutions which converge 
around the de Sitter horizon. 

(b-1) Solutions which are specified by the characteristic exponents 

Firstly, we consider the solution with the characteristic exponent v. For this, we factor 
out the singularity at z — z^ by the transformation 

/ y \ Az / ~ \ 2s + l 

^koc };s W = ^-l) Al I"- 1-— g\ ZrM , s {z\ (3-17) 

where Ai (i = 1,2,3) takes either A i+ or Then, g{ ZrjCC ,y s {z) should satisfy the 

following Heun's equation 



dz 2 + 



5 7 e 



d a + o-_z — v — a+a^ ] „ . . . 

Now we make the variable change ( = z r /z and the transformation gr Zrj00 y s (z) = ( a+ hi. s ((), 
we find a Heun's equation 

d 



d( 2 + 



a+ — a- + 1 e 7 
— 1 1 ' — 



d( 



+ 7 - "K- ! - ^± a+1 I hUO = o , (3.19) 



C(C-i)(C-«r) 

where we used the Heun's constraint in Eq.(2.17). Then, we find a solution 



9l-{z r ,oo};s 



r(v + *--% + \) (zry + g fli/( _ )B r(n + i/-(7_ + f + |) 



r(n + i/ + (7_-| + |) 



C<J 1 IV 1 -2 r 

xF ( -n - i/ + <t+ - - - -,n + u + (T+ - - + -;<r+ - <r_ + 1; — ) , 



and the other by exchanging a + with ct_ 



(3.20) 



(3.21) 



#2 ;{z r ,oo};s 9l;{z r ,oc};s\ 

They satisfy the relations fi'~^7 1 00 }. s = fi^^ oo}. s (i=l, 2), so that a combination of them 
forms a solution with the characteristic exponent v. In fact, we find 



r(<T+-e+l) 



:-)■ 



r(a_- ff+ )r(i/-(7_ + | + |)„ 



r> + <7_-! + |) 



#l;{z r ,oo};s 



rg _ r( - + - < 7_)r(»/- - + + f + |) 



(3.22) 



gives the solution. By substituting this, we find 



RV {z r ,oo};s\ z ) 



z A Hz-l) 



1 - 



1 - 



2s+l 



~ u T(n + is-a + + % + |)r(n + i/ - cr_ + f + § ) / z 



I> + -e+ 1) 

2 



n=— oo 



T(2n + 2z/ + 2) 
iv 1 lv 1 



/ cf 1 a; 1 ^ ^ ^ 2 \ 

xFn + z/ + (7 + -- + -,n + z/ + (j_-- + -;2n + 2z/ + 2;— . 



(3.23) 



is a solution specified by the characteristic exponent za Another independent solution 
with the characteristic exponent —v — 1 is given by 



(3.24) 



(b-2) Solutions which satisfy the incoming and the outgoing boundary conditions at the 
de Sitter horizon 

The solution R\ Zri00 }- s {z) is expressed by the sum of solutions which satisfy the incom- 
ing and the outgoing boundary conditions at the de Sitter horizon. We find 



R{z r ,oo};s( Z ) 



(3.25) 
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where 



~ r(n + i/-(j + + f + |)r(n + i/-(j_ + f + f; 
x 2_ 



n=— oo 



T(n + v + cr+ - | + i)r(n + 1/ + a. - % + \) \z, 

uj 3 a; 3 z 

x F[n + u-a + + - + -,n + u-a^ + - + -;2-e;l- — 



(3.26) 



and 



(z) - z A Hz Wl Z Y 3 (l Z T +1 r(1 ~ 6) V a" 
x(-J i ^ n + v + ff+ __ + _ n + 1/ + jT ___ + _ ;c; i 



(3.27) 



(V,?) The convergence region for R v , Zr O0 y s {z), etc. 

Firstly, we consider <7£/ Zj . )00 i. s (i = 1,2) in Eqs.(3.20) and (3.21). By comparing 
Eqs.(3.20) and (3.21) with (3.1), we find that the convergence region of the series for 
9i { Zr oo}-s( z ) = -^'2) i n complex z r /z plane is inside the ellipse with foci are and 
1, where with the major axis is along the real axis extending from x r to 1 — x r , i.e., 
x r < z r /z < 1 — x r . Keeping in mind that z = 1 — x and z r = 1 — x r , we find that the 
convergence region is x < 0, if we confine x to physical values. This region corresponds 
to r > r + with r + being the outer horizon. 

Now we observe that R^ Zr ^y^z) and R~[z~L}-s^. z ) are expressed by linear combination 
of 9i ;{Zr ,ccy,s( z ) (* = M), and also R {z r ,oo};s( z ) and ^{^oo};^^) are expressed by linear 
combinations of R v . ^\.Az) and R7"~^,\ Az). Thus, we conclude that the convergence 
region of all these functions are the same region, i.e., r > r + . 

(C) The matching of two solutions 

We consider the matching of the solution which is convergent for r < r' + and the 
solution which is convergent for r > r+. We make the matching in the region r + < r < r' + 
where both series are convergent. The matching can be made between two solutions with 
the same characteristic exponent. Explicitly, we require that 

R{0,1 } ; S ( Z ) = K » R \zr,00}-M- ( 3 - 28 ) 
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for the region r + < r < r' + . The proportionality constant K v is determined by comparing 
the coefficients of z p+u ~% + ^. The result is 



K, 



z p r +v * + *r( 7 )r(<7+-e + i) 



" r(p + v - % + |)r(p + + + |)r(p + 1/ + <r+ - f + ±)r(p + v + <r_ - % + ±; 



x 



X 



2 1 2' 

) n "*T(n + i/ - f + §)r(n + i/ + 5 - f + \)V{p + n + 2u + 1) 



55 r(n + z/ + f + I)r(n + i/ + 7 -f + |)(n-p)! 

p 

E a 



r(n + v - a + + | + §)r(n + i/ - cr_ + | + |) 



r(n + i/ + o+ - f + |)r(n + i/ + cr_ - f + |)r(p + n + 2i/ + 2)(p - n)!_ 

(3.29) 

The constant includes an integer p, but is independent of it. 

(D) Solutions valid in the entire region 

The solution which satisfies the incoming boundary condition at the outer horizon of 
black hole and convergent in the entire region is expressed by 



= A 



K v {s)R\ Zr ^ hs {z) + JT_ v _ 1 ( a )i2g£ } .,(z) , (3.30) 



where A s is the normalization constant. The first can be used for r < r' + with r' + being 
the de Sitter horizon and the second can be used for r > r + with r + being the outer 
horizon. 

The other solution which satisfies the other boundary condition is obtained from a 
linear combination of two independent solutions, R" n . s and 

Kut;s = (K S RL,- S y ■ (3.3i) 

which satisfies the outgoing boundary condition at the outer horizon. 

4 The Teukolsky-Starobinsky identities 

We discuss the Teukolsky-Starobinsky (T-S) identities for perturbed fields in Kerr-Newman- 
de Sitter geometries, based on Teukolsky equation defined in Eq.(2.3), although this equa- 
tion is not applied for electromagnetic field and gravitational field because they couple 
each other in this geometry. Nevertheless, we make the analysis because we want to see 
the mathematical structure of Teukolsky equation in detail and also we want to treat the 
T-S identities in the unified manner. Of course, the following analysis is valid for spin 
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1 fields which do not couple to electric charge and needless to say that it is valid for all 
fields for Kerr-de Sitter geometries. The T-S identities are expressed by 

A^) 2s A^R s = C*R_ S (T-S identity (A)), 

(V Q ) 2s R„ s = C S R S (T-S identity (B)) , (4.1) 

where C s are Starobinsky constants and Vq and Vq are defined by 



A r V 1 + a 

The T-S identities are interesting because of their mathematical structure as well as 
to fix the relative normalization between the solutions specified by s and — s. 

(A) Differential operators 

Since the solutions are expressed in terms of x or z, we rewrite the differential operator 
v q as 

x-^( x)°~ ' ' 

■ • - -i r 



dx V 1 — x 



r +- r ' -,- ( ,-i)-fi-^r a3 fi- 



(r+ - r_)(r_ - r'_) ' ' \ z r J V z^ 
x^-^^-l) 01 (1--) • (4.3) 

az \ Z r J 

The is derived by taking the complex conjugation of Vq. 

(B) Parametrizations of solutions 

Solutions are specified by the choice of A* in Eq.(2.11) so that we have to specify these 
parameters to fix them. 

(b-1) A specific choice of Ai 

Here we consider the solution which satisfies the incoming boundary condition at the 
outer horizon and thus we choose A\ = Ax-. For others, we choose A 2 = A 2+ and 
A 3 = A 3 _. Then we have 

A\ = —ax - s, A 2 = a 2 , A 3 = -a 3 - s , (4.4) 
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so that we find 



er + = 2a 2 — s + 1, <7_ = — 2ai — 2a 3 + 1, 
7 = —2a! — s + 1, o" = 2a 2 + s + 1, e = — 2a 3 — s + 1 , 
c<j = 7 + 0" — l = o" + + o"_— e = —2a! + 2a 2 + 1 . 



(4.5) 



Then, from Eqs.(3.1) and (3.30) the solution is given by the following two expressions 
which have different regions of convergence, 



A s (-x)- s - ai (l-x) a2 



x — x r 



x — x c 



2s+l 



. 1 X'p J \ 1 *^oo - 

oo 

x a n( s )F(~ n ~ v ~ fl i + a 2i n + v ~~ °i + fl 2 + 1; ~2ai - s + 1; x) , 

n=— oo 

(4.6) 

which is convergent for r < r' + (the expression (A)) and 



(4.7) 



where 



2s+l 



v F(n + z/ + a 3 - a 4 + l)r(w + z/ - Oi - a 2 + s + 1) ,^y+v+ai-« 2 
X n f^ an(Sj T(2a 2 + 2a 3 + l)r(2n + 2z/ + 2) U 



x F ( n + v + ai + a 2 - s + 1, ra + v — a 3 + a 4 + 1; 2n + 2v + 2; — ) , 



(4.8) 



which is convergent for r > r + (the expression (B)). 



(b-2) Another choice of Ai 

It may be interesting to ask the relation between solutions with different choices of 
Ai. To answer this question, we consider a generic Heun's equation, 

{ d 2 [7 5 e d cr+cr-Z — q 

) 1_ _L _i 1 1 

I dz 2 



z z — 1 z — an 



f(z)=0, 



(4.9) 



dz z(z — l)(z — an) 

with o" + +o"_— e = 7+5—1. Here we express a particular solution as Hf(an, Q] & + , <7_, 7, 5; 2) 
which is regular at z = and is normalized to be one at z — 0. By setting 

1 (4.10) 



/(*) = 1 
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the above Heun's equation is transformed to another form as 



dz 2 



7 5 2-e 

+ r + 



d 

dz 



z z — 1 z — an 

(7 + 5 - cr + )(7 + 5 - a_)z — q — 7(1 - e) 



+ 



z(z - l)(z - a H ) 



9{z) = , 



(4.11) 



where we used 1 — e = 7 + 5 — cr + — <r_. The solution which is normalized and regular 
at z = is given by Hf(an, q + 7(1 — e); 7 + S — cr + , 7 + 5 — <r_, 7, 5; z). Now we have 
two different expressions of the solution of Eq.(4.9). Both are regular at z = and are 
normalized to be one at z = 0, so that two solutions should be equal each other, 

Hf(a,q;a + ,a_,>y,5; z) 
z \ 1-6 

1 if /(a H ,g + 7(1 -e);7 + £-o+,7 + <J-<7_,7,£; 2) . 

Off/ 

(4.12) 

We apply this fact to the incoming solution at outer horizon in Eq.(4.6), which is 
expressed by 

/>"■ />"■ \ s C13 /* >■ />> \ 2s - ) - ! 



V n=— 00 / N ' 

xHf(x r , -v; <t+, o-_, 7, 5; x), 



1 - x r 



(4.13) 



where 



if/(a; r , -u;(T + ,<T_,7,<J;a;) 



E 



™ 00 1 a; 1 . 

E a„(s)F(-n- z/+ - - -,n + u + - + -;7;x) , 



(4.14) 



By using the relation between Heun's functions in Eq.(4.12), this solution is at the same 
time expressed by 



in;s 



E <(«) 



-a— 01 1 



1 — X 



where 



xHf(a H , -v' + 7(1 - e); 7 + 5 - <r+, 7 + 5 - <r_, 7, 5; x) , 



Hf(a H , -v + 7(1 - e);7 + 5 - o- + ,7 + 5 - <7_,7, 



(4.15) 



00 

E TO 

.n=— 00 



s ^/ a; 1 a; 1 

E K(s)F(-n-v+---,n + v + - + -; T ,x) 



(4.16) 
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Here the coefficients are defined by 



<(s) 



r(i/ + f-(7 + + |)r(i/ + f-(7_ + |) 



X 



r(n + 1/ + % - cr+ + §)r(n + 1/ + 1 - (t_ + §; 



(4.17) 



r(n + i/ + ct+ - f + i)r(n + v + <r_ - f + §) ' 

The difference of coefficients are due to the change of parameters in the recurrence relation 
which coefficients satisfy. This solution in Eq.(4.16) is the one which corresponds to 
A 1 = A 2 = A 2+ and A 3 = A 3+ . 

(C) Various relations 

The characteristic exponent v, coefficients a^, b v n and the proportionality constant K v 
satisfy the following relations, which are important to show that the solutions satisfy the 
T-S identities. 

(c-1) The relation between v{s) and v{—s) 
As we prove in Appendix, we have 



v(s) = v(—s) , v(s) = v(s)* . 



(4.18) 



These relation is vital for the T-S identities, because the T-S identities are differential 
transformation from the solution specified by s to that by — s. 

(c-2) The relation between coefficients and K v 

By examining the recurrence relation in Eq.(3.2), we find 



-«) = 



r(is + ai + a 2 - s + 1) 



r(n + v + a-y + a 2 + s + 1) 



T(n + v + ai + a 2 - s + 1) 



r(z/ + ai + a 2 + s + l) 
where we chose ap(— s) = ag(s). Also, by using this relation, we find 

KX-s) = v n {s) . 



a» n {s) , (4.19) 



(4.20) 



and 



K v {s) K-v-^s) r(-2ai-s + l) 



K v (-s) K.^i-s) r(-2ai + s + l) 

(c-3) Useful mathematical formula 

The following mathematical relations are needed: 



(4.21) 



l-x n 



d 

dx 



I- x r 



\ dx , 



1- X r . 



fc-1 
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dx 



(l-x) A+B ~ u F(A,B;C;x) 

T(C - A + k)T{C -B + fc)r(C) 
T(C - A)Y{C - B)T(C + k) 



(1-x) 



A+B-C-k 



F(A,B;C + k;x) , 



z A ^F(A, B; C; z) = H^1^F(A + k, B; C; z) . 



(4.22) 



(D) Teukolsky-Starobinsky identity- 
First we discuss the solution given by the expression (A) in Eq.(4.6) and show analyt- 
ically that this solution satisfies the T-S identity (A) when we take the relative normal- 
ization between R v in . s and R" n ._ s properly. In other words, we can fix the normalization 
factor for FL" n . s1 A s for s > by the T-S identity once we fix A_ s = 1. First we apply the 
operator A^.(Vq) 2s A s r to the expression (A): 



-■A, 



a 



- — (r+ - r_)(r+ - r_)(r_ - r'_) 



2s 



(-x) s ~ ai (l-x) 



S — 0,2 



s-a 3 



X 



1-X 



co 



-4s 



1 — X 



d 

dx 



~1 2s 



X — X 



1 — X 



2s+l 

x \ /-, _.\s+2a 2 



(I-*)' 



CO 



x a^(s)F(— n — v — a x + a 2 , n + z/ — a x + a 2 + 1; — 2a x — s + 1; x) 

n=—oo 
(X 

(r+ — r_)(r+ — r_)(r_ — r!_) 



2s 



X — X 



,s— ai 



(1-X) 



s— a2 



X — X r 



S — 0,3 



X — X r 



2s 



2s+l / J 



1 X y J \ 1 *^ oo 

CO 

x a p n (s)F(-n - v - a x + a 2 , n + ^ - ai + a 2 + 1; -2a x - s + 1; x) 



71= — OO 

(r+ — r_)(r+ — r_)(r_ — r'_) 



0/ 



2s 



X 



r(-2 ai -g + i) 

r(-2ai + s + l) 



r(i/ + ai + a 2 + s + 1) 



- rt in;— s' 



T(z/ + ai + a 2 - s + 1) 
where we used relations in Eqs.(4.22) and (4.19). Thus A s is fixed by 

„2 n2s 



(4.23) 



A s 



c: 



a(r+ — r_)(r+ — r_)(r!_ — r_) 



x 



r(-2ai + s + i) 
r(-2ai - s + 1) 



r(z/ + ai + a 2 - s + 1) 



r(z/ + a! + a 2 + s + l) 
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(s > 0) , 



(4.24) 



provided A_ s = 1 (s > 0). 

Next we show that the expression (B) of the solution in Eq.(4.7) satisfies the T-S 
identity (A). To the end we observe for s > 



a 



— (r+ - rJ)(r' + - r_)(r_ - r'_) z s - a \z - 1) 



a 
x ( ! 



2s 



~ \ s— 0,3 



-4s 



1 - 



Z, 



dz 



2n 



1 - 



z 



oo 



-2s+l 



Z 



S + 2CL2 



^ u T(n + v + a 3 - a 4 + l)T(n + z/ - a x - a 2 + s + 1) 
X n f=' oo a "^ r(2a 2 + 2a 3 + l)r(2n + 2z/ + 2) 

n+!/+ai— a2 



X 



a 



— F [n + v + ai + a 2 - s + 1, ra + i/ - a 3 + a 4 + 1; 2n + 2v + 2; — 



— (r + — r_)(r\_ — r_)(r_ — r'_ 



2s 



xz s - a2 ( z - iy- ai i — 



S— 03 



1 



-2s+l 



2s 



I d V" s+2a2 \^ v I \ T( n + v + a 3 - «4 + l)T(n + v - a\ - a 2 + s + 1) 
X ~ ' Z ni^oo r(2o 2 + 2a 3 + l)r(2n + 2z/ + 2) 

x ( — | F I n + zy + ai + a 2 - s + 1, ra + f - a 3 + a 4 + 1; 2n + 2v + 2; — 



dz 

z 



Z r 



a. 



— (r + — r_)(r' + — r_)(r_ — r'_) 



2s 



T(v + ai + a 2 + s + 1) 



r(z/ + ai + a 2 - s + 1) 



T3V 

- rt {z r .,oo};-s ) 



(4.25) 



where we used relations in Eqs.(4.22) and (4.19). 

By exchanging v with —v — 1 and using the fact that 2s is integer, we find from 
Eq.(4.25) that 



a 



— (r + — r_)(r' + — r_)(r_ — r'_) 



2.s 



r(i/ + a! + a 2 + s + 1) 



Y{y + ai + a 2 - s + 1) 



X1 {2 r .,oo};-s 



(4.26) 



Thus we find 



A;(^) 2s AX ;s = is £(r + -r_)(r' + -r_)(r_-r'_) 



2s 



X 



r(-2ai -s + 1) 
r(-2ai + s + 1) 



r(z/ + ai + a 2 + s + 1) 



r(z/ + ai + a 2 - s + 1) 
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x K v {-s)R\ ZrM; _ s {z) +K_ v _ 1 {-s)R- { l-^_ s {z) 
= CI [K v {-s)R\ Zr ^_ s {z) + K. v . 1 (- 8 )^J oy ._ a (z) 



(4.27) 



where we used the relation in Eqs.(4.21). 

Next we examine the T-S identity (B) by using another parametrization of the solution 
given in Eq.(4.15). By applying the differential operator (T>q) 2s to this expression, we 
find 

/ \ s—2a* r oo n r oo i _J - r „ ^> 1 2s 

\2s rtv 



f -x r y-2°s 

\1 — X r J 


oo 




oo 


-1 - 


E <(-*) 

_n=— oo 




E K(s) 

_n=—oo 





(r+ — r_)(r_ — r^_) 



x(-x) ai (l -x) 



«2 



a 3 



1 J dec 



d 



~1 2s 



1 — X 



2s+l 



x b u n (-s)F(-n - v - ai + a 2 ,n + ^ - ai + a 2 + 1; — 2ai + s + 1; x) 



/ -X r \ 2s 
\l-X r J 


oo 




oo 


-1 r 


E <(-*) 

_n=— oo 




E <(«) 

_n=— oo 





r + — r_ 



(r+ — r_)(r_ — r'_ 



2.s 



r(-2a 1 + g + l) 



(4.28) 



r(-2ai -s + 1 

where we used relations in Eqs.(4.22) and (4.20). 

Thus the T-S identity implies that the following identity among the sums of coefficients 
should be satisfied, 



E <(-s) = \c s 



~1 2s 



a(r' + — r + )(r_ — r'_) 



T(u + ai + a 2 - s + 1) 



r(z/ + ai + a 2 + s + 1) 



E <(*) • 

i=— oo 

(4.29) 



5 Summary and discussions 



In our previous paper [8] , we showed that both angular and radial equations for the Teukol- 
sky equation in Kerr-de Sitter (for all massless fields) and Kerr-Newman-de Sitter (mass- 
less fields except electromagnetic and gravitational fields) geometries are transformed to 
Heun's equation and derived analytic solutions in the form of series of hypergeometric 
functions. The solution of the radial equation presented in Ref.8 is convergent for r < r' + 
with r' + being the de Sitter horizon. In this paper, we constructed the analytic solution 
of the radial equation which is convergent for r + < r with r + being the outer horizon. By 
matching these solutions in the region where both solutions converge, i.e., r + < r < r' + , 
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we obtained the analytic solution valid for the entire region of r. We showed analytically 
that our solution satisfies the Teukolsky-Starobinsky identities which are quite nontrivial 
differential transformations. 

Although the use of our solution in this paper is less clear than that for the Kerr 
geometry case, we believe that our solution will be relevant to the physical world, espe- 
cially in the early stage of the universe where the cosmological constant plays important 
role. By using our solution, we can derive the analytic expression of the decay rate for 
emission of massless particles from the Kerr-de Sitter or the Kerr-Newman-de Sitter black 
hole [11], which can be used to examine a correspondence between quantum gravity on 
anti-de Sitter space and a conformal field theory defined on its boundary [12,13], 
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Appendix 

(a) The proof that A~ S R*_ S is a solution 

Form the radial Teukolsky equation for R s , we derive the equation for A s r R s . We find 

eQr \ 2 ( ' eQr \ dA r 



a ^ d a i 9 d 1 
A — A 1 

r dr r dr- A r 



(l + a) 2 (K-^) -i8(l + a )(K- 



1 + a J dr 



2a 



+ 4is(l + a)uor - — (s - l)(2s - l)r 2 - s(l - a) - 2iseQ - \ s }A s r R s (r) = 0(A.l) 



By using the relation X s = A_ s , we find that this is the equation which (R_ s )* satisfies. 
By changing s to — s, we conclude that (A~ s i?_ s )* satisfies the same Teukolsky equation 
as R s does. 

(b) The proof of v(s) = v(—s) 

With this choice of Aj in Eq.(4.4), coefficients in the recurrence relations are expressed 

by 

(n + v + a x - a 2 + l)(n + v + a 3 - a 4 + l)|n + v + a x + a 2 + s + 1| 2 



<(s) 



2(n + z/ + l)(2n + 2z/ + 3) 

(gi - a 2 )(a3 ~ a 4 )[s 2 + (gi + a 2 )(a3 + a 4 )] /I _ \ , w n 
2(n + i/)(n + i/ + l) +l v2 ^r^n + ^Kn + i/ + ij 



+ 



/ \2 



s 2 - aia 3 + a 2 a 4 + - + x r (a\ + a 2 ,) 



2 (r + — r_)(r+ — r'_ 
2a 4 (l + a) 2 (r_ — r+) 2 (r_ — r'_) 2 (r + — r'_)(r+ — ri.) 



a 2 P 



r + — r. 



x | — uj r + (r + r_ — 2r_r' + + r + r' + ) + 2auj(au — m)r + (r_r' + — r + ) 

eQ 

—a 2 (au — m) 2 (2r + — r_ — r+) + -^—^ — [ cjr + 2 (r + r_ + r 2 — 3r_r^ + r+r^) 
— a(acj — m)(r + r_ — 3r 2 + r_r' + + r + r^) ] + ^ r +( — r + + r - r +) | 



2r 



/ 2 



(r + — r_)(r+ — r'_) 



+ Xr 



a(r + — r_)(r+ — r'_) 



(n + v — ai + a 2 ) (n + v — a 3 + a 4 ) |n + z/ + a\ + a 2 — s\ 
2{n + u){2n + 2u - 1) 



(A.2) 



And we see that 



<-i(*)tM = 



\n + u + a ± — a 2 \ 2 \n + u + a 3 — a 4 \ 2 \(n + u + a ± + a 2 ) 2 — ° 2 ' 2 
4(n + z/) 2 (2n + 2i/ - l)(2n + 2v + 1) 



-, (A.3) 
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Now we observe that the transcendental equation which determine v given in Eq.(3.6) 
contains coefficients as /3^(s) and a special combination ctn-i( s )ln( s ) ■ From expressions 
in Eqs.(A.2) and (A. 3), we find that both (3^(s) and «^_i(s)7^(s) are even functions of 
s. Furthermore, we see that the transcendental equation for v is an equation of real 
coefficients. Therefore if v(s) is a solution, then v{—s) is a solution, and also v(s)* is a 
solution. By the uniqueness of the solution, we conclude that the solution v(s) is an even 
function and a real function of s. 
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